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I. INTRODUCTION 


Pulsed ultrasonic techniques are standard methods for 
detecting and characterising defects in materials and 
devices .*” 4 To accurately determine the defect parameters 
and structural integrity of a multilayered structure using 
ultrasonic techniques, the evaluator must understand the 
wave propagation in the structure and identify the 
appropriate signal of interest. Interpreting the signal 
response is relatively easy with a single-layered specimen, 
but becomes complicated as number of layers increases, 
because of multiple pulse reverberations. 

We have derived generalized transit time and pulse 
amplitude formulas for modelling the pulse reverberations in 
multilayered structures. The equations can be programmed 
into computer routines and using the appropriate parameters, 
will calculate the time and echographic amplitude response, 
which can then be used to identify the critical structural 
interfaces . 

The equations are applicable for both the reflected 5 
and the through-transmitted modes. The reflected mode is 
not applicable in materials and structures that are highly 
attenuating because the reflected signals are either buried 
in the noise or too weak to be detected. As long as the 
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material is accessible from the side opposite the 
transmitter, through-transmission is a viable alternative 
for evaluating lossy materials. 

In this paper, we show how the generalized formulas for 
through-transmission time and pulse amplitude response were 
derived, and present the results of laboratory experiments 
to verify them. Finally, we discuss the validity of the 
equations and their applicability to ultrasonic 
nondestructive evaluation. 

Acknowledgement. The first author was previously with 
Martin Marietta Laboratories. This work was initially 
supported by the Martin Marietta Laboratories' (MML) 

Directors Discretionary Fund Program and partially supported 
by the MML Independent Research and Development Program. 


II. THEORETICAL DERIVATION AND ANALYSIS 

Stress-wave propagation in multilayered structures is a 
very complex phenomenon. In addition to complications 
caused by the normal modes of propagation, other factors 
such as surface finish, interface quality, diffraction, 
scattering, attenuation, focusing, and mode conversion also 
strongly influence the transmitted and reflected signals. 
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Because these factors are mostly material and process 
dependent , they are very difficult to model mathematically. 

In the derivation below, we consider only the pure mode 
through- transmission of a pulse that is normally incident to 
a n-layered structure, which is bounded by two semi-infinite 
media on each end. Figure 1 illustrates wave propagation in 
a typical n-layered structure (d is the thickness of the 
layer and a is the attenuation factor, with the subscript i 
for the ith layer) . The acoustic impedance is given by Z » 
pv, where p is the density of the medium and V is the sound 
velocity. We assume that all adjacent layers have mismatched 
acoustic impedances, i.e., Z^ ^ + i» 

The reflection and transmission coefficients for an 
ultrasonic pulse propagated from medium i to medium j are 
given by R^j * (Z j-Z^) / (Z j+Z^) and T^j « (2Z j) / (Z j+Z^) , 
respectively. For a single-layered structure, the analysis 
is straightforward. The time of arrival for the nth pulse 
reverberation can be expressed as 


t (n) = (2n-l) t j * (2n-l)d l /v i 


( 1 ) 
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where t^ * d^/v^ i® the time it takes the pulse to travel 
through thickness, d lf with the wave velocity of v x . The 
pulse amplitude associated with the specific pulse that 
arrives at t can be written as 


n-l n-1 - (2n-l) d.a, 

A (n) = A 0 T 12 R 10 R 12 C 


( 2 ) 


where Aq is the initial pulse amplitude. As the number of 
layers in a structure increases, the echo response becomes 
more complex because the number of reverberations increases 
and the echo become superimposed. Each transmitted echo 
E (k^, ... , k , ... , kj^) is indexed by the number of txmes rt 
reverberates in each layer, and a specific echo is 
represented by a specific set of indices, k^ The indices 
are used to determine the number of echo elements that form 
the echo. Thus, a pulse response at a given moment may be 
composed of many echoes that have the same time of arrival. 
This phenomenon is demonstrated in Figure 2. 

Depending on the phase differences of these echoes, 
there is either constructive or destructive interference. 
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which causes an irregular increase or decrease in amplitude, 
and which produces local minima and maxima. The formula to 
compute the number of echo elements is based on the binomial 
coefficient that determines the number of outcomes that can 
occur within a given set. A multiplication factor based on 
Pascal's triangle (which can also be derived based on 
binomial coefficients) specifies the sire of a subset of 
echo elements that follow similar reflection and 
transmission patterns through the material. The formula is 
given by — 


N = 


n 

j=i 


nnn(k r k w ) 

I 


i=l 


k. -1 
■ 


j-1 



m! 


n! (m-n)! 


(3) 


Based on the above analysis, we can derive generalized time 
and amplitude equations for the through-transmission mode. 
The time response is derived via a linear summation of the 
time required to travel through each layer. The 
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mathematical expression for the transit time response of a 
pulse can be written as — 


t 


^rV 



X <*■-•> ‘i 


(4) 


where k^ * l f 2 , . . . n, and t^ - d^/v^ i® the time required 
for the pulse to travel through the ith layer. 

The corresponding amplitude at each calculated time is 
computed by observing the possible paths the pulse might 
travel in the given amount of time and is expressed using 
the k^ indices from which the echo elements are determined. 
Each echo element follows a definite pattern of reflections 
and transmissions at the material interfaces before it exits 
to the receiver. By integrating these reflection and 
transmission coefficients with the results of Equation (3), 
the echo response at a given time is defined for a given set 
of indices. The equation for the pulse response amplitude 
can be expressed as 
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(k T k » k n * 


( 5 ) 


n»n+l 


k, • 1 ' 


n T .,-, 

1*2 


10 


R . e 

n,n+l 


- I (2k- -1) d ,a 
i«l i 


MV, 


n * 

i=l 




* n 



Aq is very difficult to measure experimentally because 
it depends on several factors, including the transducer 
characteristics, pulse driver, and coupling efficiency. 
However, we can normalize the expression for the n-layered 
structure to kj ■ 1 and k n ■= 1, i.e., (n-2) layers. Thus, 

the expression for Equations (4) and (5) reduced to -- 


n- 1 


n - 1 


•u.k-.-k r i>- *,+ E ( 2 k r !) V. - *. <«> 

2 H' 1 i=2 i=2 
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and 


0-^ W> 


( A 0 T 12 T n,n+1 


- d ,a , - d a 
e 11 " ■) * 


n-l 


T n-l,n ‘-** 1 

R R e 

n-l 21 n-l,n 

n t u ., 

i-3 


M*i + l 

n r 

A A M+l 
i=2 


n-2 

n 

i*2 


j-i \ m / \ j-i 


* 

J (7) 

j 


R 2 

K i.i+1 


Equations (6) and (7) are the same expressions as the 
formulas specifically derived for the (n-2) —layered 
structure. Therefore, the first through-transmitted pulse, 
which can be measured experimentally can be used as the 
normalization factor. 

The time and amplitude data can be plotted to provide a 
pictorial presentation of the effect of multiple 
reverberations that can be used to predict the contribution 
of significant echoes from the interfaces between layers. 
Based on the plot, an analysis can be made to determine 
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which transmitted pulses are of interest for a specific 
inspection, and signal processing algorithms can be 
developed to eliminate unwanted reverberations . g 


III. LABORATORY EXPERIMENTS AND RESULTS 


Ultrasonic immersion tests were performed to verify the 
validity of Equations (4) and (5) for one- and three-layered 
structures. A benefit of using immersion tests is that it 
eliminates the error introduced by the thin coupling 
required in the contact tests. As described in Equations 
(€) and (7) , we first normalized the measurements for a 
three-layered structure to t lm ^ and A^ m ^ and reduced it to 
verify the equations for a one-layered structure. 

Similarly, we normalized the measurements of a five-layered 
structure to and A^jj^, and reduced it to three 

layers . 

The equations for the three-layered structure are as 
follows : 


(k i i = < 2k r ^ *1 + ^2* ^ l 2 + ^ k 3* ^ l 3 


(8) 
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and 


A (k ,k Jc ) * A 0 
k 12 3 1 


34 


T T 

21 32 


k -1 k -1 k +k k +k 
P 1 p 3 p 1 2 P 2 3 * 

10 34 K 12 23 * 


- Ok , -1) d, „ (21 2 d 2 „ Bk, -1) d, a , 

V V 


min(k j Jc 2 ) 

l <- 


V‘ 


1) 


inin(k 2 *k 3 ) 

Y 

" (-D 


Vi 


j=l 


i-1 


k,-l \ / 


j-1 /\ 



(9) 


The reduced expressions for Equations. (8) and (9), i.e., 
for one layer, are 


*<1 . k 2 .l ) 


+ (2k 2 - 1) t 2 + t 3 


( + t 3 ) +(2k 2 - l)t 2 (10) 


and 
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0 . k 2 .l ) 


-da -da 

11 i 3 

< A 0 T 34 T 12 e > 


k 2-‘ 

T 23 R 2? 


k - 1 
*23 * 


(2k 2 -l)d 2 o 2 


The time-domain offset (t^ + tg) and normalization amplitude 
A(j.ii) can be measured experimentally. 

Three experiments were performed to verify the time and 
amplitude formulas: two experiments using a three-layered 
structure (water/0 . 102-cm glass/water, water/0 . 131-cm 
glass/water) and one using a five-layered structure 
(water/0 . 318-cm aluminum/O . 102-cm water/0. 620-cm 
aluminum/water) . The experimental setup is shown in Figure 
3. The measurements were taken with a Krautkramer-Branson 
USIP 12 pulser /receiver and a Tektronix 7904 oscilloscope. 
Two 10 -MHz transducers were used to transmit and receive the 
ultrasonic pulses. The specimens were immersed in a water 
bath approximately 1.27 cm from the transducers. 

The acoustic properties of the materials used in these 
experiments are listed in Table 1 . 


11 



Table 1 - Acoustic properties of the materials used in 


experiments 


Material 

Glass 

Aluminum 

Water 


Acoustic Impedance 
(kg/m 2 /s) 

14500 

17060 

1483 


Sound Velocity 
(m/s) 

5710 

6320 

1483 


Figures 4 and 5 show the oscilloscope traces of the 
pulse reverberations for a single-layered (0.102-cm glass 
plate) sample and a three-layered (0.318-cm aluminum/0 . 102- 
cm water/0. 620-cm aluminum) composite sample, respectively. 
The theoretical and laboratory data for these three cases 
are compared in Figures 6-8. The theoretical calculations 
of transit time responses predicted using Equation (4) agree 
well with the experimental measurements. The theoretical 
calculations of the echo amplitude response [Equation (5) ] 
also agree well with the experimental results for the 
single-layered samples [Equations (10) and (11) , and Figures 
6 and 7], but there are discrepancies between the 
theoretical and experimental results for the multilayered 
samples. We attribute some of the error 1) to the limited 
accuracy of the scope, 2) to imperfections in the surface 
finishes, 3) to beam spread and dispersion, 4) to imperfect 
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alignment of the sensor, and 5) to the superimposing of 
multiple pulses. Although the experimental measurements are 
far from ideal, they are adequate to verify the general 
principles of pulse response as predicted by Equations (4) 

and (5) . 


IV. CONCLUSION 

It is evident from the results that the generalized 
equations for through-transmission pulse response in 
multilayered structures, which are derived here for the 
first time, are indeed valid. We have verified them using 
two approaches. First, we reduced the equations to specific 
forms of single- and three-layered structures and found that 
the resultant formulas are of the same form as those derived 
for each cases individually. We also verified the equations 
experimentally, with generally good agreement between the 
predicted and experimental measurements of transit time and 
amplitude response. 

These equations are invaluable tools for modeling 
ultrasonic pulse responses and for practical inspections. 
They are applicable to both thick and thin specimens, which 
satisfies the condition d i »wavelength, and can be applied 
to small electronic components as well as to large 
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structures. In addition, they can be easily incorporated 
into an expert NDE system for evaluating the interfaces of 
multilayered structures . 
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Figure 1. Wave propagation in a typical n-layered 
structure. 
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For n - 2 , 




Figure 2. Different wave passages with the same pathlength 
and arrival time. 
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Figure 3. Equipment setup for laboratory experiments. 
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Figure 4 . 


Oscilloscope trace of the pulse reverberation 
for a single 0 . 102-cm-thick layer of glass. 
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Figure 5. Oscilloscope trace of the pulse reverberation 

for a three -layered (0.318-cm aluminum/ 0 . 102-cm 
water/0 . 620 -cm aluminum) sample. 
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Figure 6. Comparison of the theoretical prediction 

[Equations. (10) and (11)] of pulse amplitude 
and arrival time with experimental results for 
single 0 . 102-cm-thick layer of glass. 
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Figure 7, Comparison of the theoretical prediction of 
pulse amplitude and arrival time [Equations . 
(10) and (11) ] with experimental results for a 
single 0 . 131-cm-thick layer of glass. 


22 




TIME (^s) 


Figure 8. Comparison of the theoretical prediction of 

pulse amplitude and arrival time [Equations. (4) 
and (5) ] with experimental results for a three- 
layered aluminum/water/aluminum sample. 
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